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Abstract 

The problem of resonant generation of nonground-state condensates is addressed 
aiming at resolving the seeming paradox that arises when one resorts to the adia- 
batic representation. In this picture, the eigenvalues and eigenfunctions of a time- 
dependent Gross-Pitaevskii Hamiltonian are also functions of time. Since the level 
energies vary in time, no definite transition frequency can be introduced. Hence no 
external modulation with a fixed frequency can be made resonant. Thus, the reso- 
nant generation of adiabatic coherent modes is impossible. However, this paradox 
occurs only in the frame of the adiabatic picture. It is shown that no paradox exists 
in the properly formulated diabatic representation. The resonant generation of di- 
abatic coherent modes is a well defined phenomenon. As an example, the equations 
are derived, describing the generation of diabatic coherent modes by the combined 
resonant modulation of the trapping potential and atomic scattering length. 
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1 Origin of adiabatic paradox 



Nonequilibrium ultracold Bose-condensed gases are well described by the time-dependent 
Gross-Pitaevskii equation [1-10], which is a nonlinear Schrodinger equation 

^m = H[^|;,t]m■ (1) 

Here ?i = 1 and, to simplify the notation, the condensate wave function '0(t) is assumed to 
be a vector in spatial variables and the nonlinear Hamiltonian i^f [-0, t] is a matrix in these 
variables. The ovcrdot means the differentiation with respect to time t. The temperature 
is close to zero, so that all atoms are supposed to be in the coherent Bose-condensed 
state. Generally, the condensate wave function is normalized to the total number of 
condensed atoms. If this function is denoted as rj{t) , it is always possible to introduce, by 
means of the relation 

the function i/j{t) that is normalized to one, 

mm' ^<mm)>- i- (2) 

It is the latter function that is assumed in Eq. (1). The system of trapped atoms is 
subject to the action of an externally induced modulation field, so that the nonlinear 
Hamiltonian H[ilj,t\ depends on time directly as well as through the function ip{t). 
It is possible to consider the adiabatic eigenvalue problem 

H['i/;^,t]Mt)^E^{t)Mt) , (3) 

treating time as a fixed parameter. Then, of course, the eigenfunctions ipnit) and the 
eigenvalues En{t) depend on time. The multi-index n labelles the eigenfunctions and 
eigenvalues. The energy spectrum, for trapped atoms is discrete. The eigenfunctions can 
be normalized to one: 

\\Mt)\\' =<Mt)\Mt)>- 1- (4) 

These functions correspond to the adiabatic coherent modes. 

The nonlinear Eq. (1) can possess different solutions. A wide and important class of 
solutions can be represented as the expansion over the adiabatic coherent modes: 

i;{t) = ^a„(t) ex.p{ix„{t)} ipn{t) , (5) 

n 

where the phase 

Xn{t) = 6n{t) + Cn{t) (6) 

is the sum of the dynamic phase 

5n{t) = - f E^{t')dt' (7) 

and of the geometric phase 

Cn{t)=l f <ijn{t')\'^n{t')>dt' . (8) 
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Formula (5) does not exhaust all admissible solutions of the nonlinear Schrodinger equa- 
tion (1). But here we are interested in the class of solutions that are representable as 
expansion (5). 

The problem of dealing with the nonlinear Hamiltonian is that the form if['0n) t] 

is not a Hermitian operator since 

< i^mit) I H[iljn,t]Mt) > + < ayi^nAi'mit) \ > ■ (9) 

Therefore, the eigenfunctions ■0m (^) a-nd for m ^ n, generally, are not orthogonal, 

^mn ■ 

These functions do not necessarily compose a complete basis. However, in the space that 
is a closed linear envelope Span{il)n{t)}, they form a total normalized basis, so that the 
functions from this space can be represented in the form of expansion (5). Substituting 
the latter into Eq. (1) yields 

^ [arn{t) < i^„{t) I Vm(^) > + am(^) < i^n{t) \ i^m{t) > (1 - Smn) + 
m 

+ ia^it) < iJnit) I H[ij, t] - E^{t) I ^l:^{t) >] exp{^x™(^)} = . (10) 

This equation is difficult to simplify because of the nonorthogonality of the basis 

If the time dependence in the Hamiltonian H[%l),t] enters through an external alter- 
nating field with a fixed frequency a;, then the latter cannot be tuned to resonance with 
any of the time-dependent transition frequencies E^{t) — En{t). Moreover, the eigenen- 
ergies enter expansion (5) not explicitly but through the dynamic phases (7). Thus, no 
resonance condition can be defined, making the resonant generation of adiabatic coherent 
modes impossible. It is this argument that one raises against the possibility of generating 
nonground-state condensates of trapped atoms. 



2 Paradox-free diabatic representation 

Now we show that no paradox arises in a properly formulated diabatic representation. 
Let the Hamiltonian be a sum 

i/[^,t] = i7oM + ^[^,t], (11) 

where the first term describes the system of cold trapped atoms, while the second term 
contains the direct time dependence from externally induced modulation fields. 

The diabatic coherent modes are defined as the solutions to the time-independent eigen- 
problem 

Ho[iPn\ fn = E^fn ■ (12) 

The definition of the diabatic coherent modes as time- independent, stationary solutions 
makes them principally different from the adiabatic modes of Sec. 1. It is this definition of 
coherent modes that was introduced in Refs. [11-13]. Particular examples of such modes 
are vortices [14], though many other types of modes are admissible [11-13]. The lowest 
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diabatic coherent mode corresponds to the usual ground-state Bose-Einstein condensate. 
While the excited diabatic coherent modes represent nonground-state condensates. 

The nonlinear Hamiltonian i/gi'/'n] is not Hcrmitian. Hence the cigcnmodcs ip^ and 
(fni with m ^ n, are not orthogonal, although all of them can be normalized to one: 

llV^nlP =<^nWn>= 1- (13) 

But again, in the space that is a closed linear envelope S'pan{(/9„}, the family forms 
a normalized total basis. In this space, the corresponding class of solutions to Eq. (1) 
can be represented as the expansion 

^{t) = E Cn(^) QM-iEnt)^n (14) 
n 

over the stationary coherent modes. 

Substituting expansion (14) into Eq. (1) results in the equation 

Cm{t) <'Pn\Vm> exp(-ia;^„i) = 

m 

= ^rn{t) < \ H[^, t] - En \ <^m > exp{-iuJmnt) , (15) 
m 

in which 

^mn = — En . (16) 

The latter expression defines transition frequencies that do no depend on time. Conse- 
quently, a given frequency cu of an external alternating field can be tuned to the resonance 
with one of these transition frequencies. 

Equation (15) can be simplified assuming that the coefficients c„(t) are slow functions 
of time as compared to the fastly oscillating exponentials, so that 

Cnit) 



E„ 



< 1. (17) 



Then, we can use the averaging method [15] and the scale separation approach [16-18]. 
Employing these techniques, we average Eq. (15) over time, defining the time averaging 
as 

{/(t)}„„ = Jim i [ fit)dt, (18) 

and treating the coefficients Cnit) as quasi-integrals of motion. The averaging of the 
exponentials gives 

•{exp( i^mrJ^}av — ^mn ■ 

Introducing the notation 

l^nm{t) = {<ipn \ H[%l),t] - Em \ (fm > exp{-iLJmnt) } av , (19) 

we come to the equation 

iCnit) = Y l^nmit) Cm{t) . (20) 
m 

The solutions c„(f) define the time variation of the fractional mode populations |c,i(t)|^. 

The principal difference of the diabatic representation, described in this section, from 
the adiabatic picture of Sec. 1, is as follows: 
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• There exist well defined stationary energy levels of the diabatic coherent modes. 

• The frequency uj of an external field can be tuned to the resonance with one of the 
transition frequencies ojjnn- 

• The diabatic coherent modes, being the solutions to the stationary eigenproblem, 
do not depend on time. 

• The equation for the coefficient functions of expansion (14) can be simplified by 
means of the averaging techniques. 

These essential points make it possible to consider well defined resonance conditions. 
Consequently, it is feasible to realize the resonance generation of the diabatic coherent 
modes. 



3 Generation of coherent modes 

To specify the above equations, let us consider the atoms interacting through the local 
potential 

$(r) = $o5(r) ($0 = 47r ^) , (21) 

where m is atomic mass and a^, scattering length. The nonlinear Hamiltonian takes the 
form ^ 

H[^A = - ^ + U{v,t) + iV<|.(t)|^(r,t)|2 , (22) 

where the term 

U{v,t) = U{v) + V{v,t) (23) 

consists of a trapping potential U{r) and of an external driving field V{r^t). The inter- 
action part 

$(t) = $o + £(t) (24) 

contains the static interaction potential $o, defined in Eq. (21), and an additional term 
describing a possible modulation of the scattering length by means of the Feshbach reso- 
nance technique [19]. 

The generation of coherent modes by modulating the trapping potential, as in Eq. 
(23), their properties, and different applications were considered in Refs. [11-13,20-35]. 
The analytical treatment, using the averaging method, was compared and found to be in 
good agreement with the direct simulation of the Gross-Pitaevskii equation [36,37]. The 
possibility of the resonant creation of coherent modes by modulating the atomic scattering 
length using Feshbach resonance was considered in Ref. [38]. This method can be pre- 
ferred for those atomic species that demonstrate high tunability of their scattering length 
[39]. In the present paper, we take into account both these ways, simultaneously modu- 
lating the trapping potential as well as the scattering length. By combining both these 
techniques can give additional advantage for generating nonground-state condensates. For 
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instance, the combination of these two techniques can be employed for generating two dif- 
ferent excited coherent modes, in addition to the ground-state one, or for enhancing the 
generation of the chosen mode. 

Hamiltonian (22) can be presented in form (11), with 

^oM = - ^ + U{r) + iV$o|^(r, t)\' (25) 

and 

V[i;,t] = Vir,t) + Ne{tmr,t)\' . (26) 
Eigenproblem (12) becomes 

HqM <^n(r) = En (p„{r) . (27) 
And expansion (14) takes the form 

i;{r,t)^J2 Cn(i)e-'^"Vn(r) . (28) 

n 

Then we should follow the consideration of Sec. 2. 

To explicitly accomplish the averaging procedure of Eq. (18), we have to concretize 
the trap modulation and the type of the interaction alternation. We take the general form 
of the oscillating trapping modulation 

V{r,t) = Vi{r) cos ujt + V2{r) sin ujt (29) 

and the similar type of the interaction alternation 

s{t) = Si cos cut + 82 sin ujt . (30) 

The oscillation frequency is taken the same for both modulations, so that the same 
nonground-state coherent mode be excited. In principle, two different frequencies could be 
taken for modulations (29) and (30), if we would wish to excite two different nonground- 
state modes simultaneously. The frequency cu is tuned to the resonance with a chosen 
transition frequency 

= E2-Ei. (31) 

If the generation procedure starts with an equilibrium Bose-Einstein condensate, then 
El is the lowest energy of atoms in the trap, corresponding to the condensate chemical 
potential [10]. The resonance condition reads as 



Au 



UJ21 



< 1 {Aio = u)- UJ21) ■ (32) 



For what follows, it is convenient to introduce the notation for the transition ampli- 
tudes 

amn = N^o j \^m{r)\^ [2|<^n(r)|^ - |^m(r)|^] dv , 



[5mn = / V*ra{v) [Vi{v) - lV2{v)] ^„(r) dv 
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7n = N{ei - 182} j ipl{v) |(/7„(r)|V2(r) dr . (33) 

The first of them is the amphtude due to atomic interactions, the second one is caused 
by the trap modulation, and the third amphtude is due to the interaction oscillation. To 
satisfy condition (17), these amplitudes are to be small, such that 



UJ-n 



« 1, 



(3n 



« 1, 



In 



OJr. 



« 1, 



where uJmn 7^ 0. 

Substituting expressions (29) and (30) into Eq. (20), under the resonance condition 
(32), yields the evolution equations for the coefficient functions = Cn{t) defining the 
fractional mode populations 

p^{t) = |c„(t)|2 . (34) 

Under the resonance condition (32), the system of equations reduces to only two equations 
for the considered coherent modes: 

■ I |2 I i2, I i2 . n \ i.Auit , 1 

t 



dt 



ai2|c2pci + ^ (27i|ci|2 + 72|C2|2 + A2) c^e''^'^' + ^ Yi^cj e"^^-* , 
a2i|cipC2 + I (272IC2I' + 7!|cip + ^;,) cie-*^-* + ^ j^clcl e*^"* . (35) 



dt ' ' 2 V ^' ' - 2 

Solving these equations gives the dynamics of the fractional mode populations (34). 

Here we have considered the case of zero temperature and asymptotically weak atomic 
interactions, when the whole atomic cloud is in a coherent state described by the Gross- 
Pitaevskii equation. It is possible to extend the consideration to the case of finite temper- 
atures and interactions, when, in addition to the fraction of condensed atoms, there exists 
a fraction of uncondcnscd atoms. This can be done by invoking, for instance, a stochastic 
variant of the Gross-Pitaevskii equation [40-42] or by employing the fully self- consistent 
theory [43-45]. The possibility of generating nonground-state condensates, even at finite 
temperatures and interactions, is due to the resonant nature of the generation procedure; 
since the energy levels of the coherent modes in a trap are discrete, while the spectrum 
of uncondensed atoms is continuous [43] . 

Equations (35) describe the dynamics of the guiding centers that correspond to the 
first approximation of the averaging techniques [15-18]. In this approximation, under the 
resonance condition (32), the system of equations for the functions c„(i) reduces to only 
two Eqs. (35), if at the initial time t — other modes were not populated, so that 

c„(0)=0 (n^l,2). 

Then c„(t) = 0, if n ^ 1,2, for all t > 0. Generally, there also exist nonrcsonant 
transitions between the coherent modes. These nonresonant transitions can be taken into 
account by the higher-order approximations of the scale separation approach [16-18]. In 
the higher-order approximations, the mode amplitudes c„(t) for the nonresonant modes, 
for which n 7^ 1,2, become nonzero, even if at the initial time these modes were not 
populated. However, the amplitudes of these nonresonant modes remain small during the 
time interval < t < tres, when they can be neglected. But for longer times t > tres, 
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the nonresonant modes cannot be neglected. Thus, the resonant generation of coherent 
modes is hmited by the time tres, which was estimated in Refs. [25,43]. By the order 
of magnitude, the resonance time tj-es can be made comparable to the lifetime of atoms 
in a trap, provided the absolute values of the transition amplitudes |q;i2|, |/9i2|, |7i|, and 
I72I are much smaller than the transition frequency CU21 ~ cu. Thus, Eqs. (35) are valid 
only for t < tres- But the time tresi being of the order of the atom lifetime in a trap, is 
sufficiently long for realizing the process of the resonant mode generation. 

In conclusion, we have demonstrated that the appearance of the adiabatic paradox, 
precluding the generation of the adiabatic modes, exists solely in the adiabatic picture. 
The adiabatic representation is unsuitable for the case of resonance. However, the diabatic 
representation, developed in Sec. 2, contains no paradoxes and is perfectly appropriate 
for describing the resonant generation of coherent modes characterizing nonground-state 
condensates. Such a generation can be realized by the resonant modulation of cither the 
trapping potential, or the scattering length, or by the combined modulation of both. 

The appearance of the resonantly generated nonground-state condensates can be ob- 
served, for instance, by means of collective light scattering [46-49]. Another way of their 
registration is through the time-of-flight experiments, as is discussed in Ref. [50]. 
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